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Spherical neutron star models are studied within tensor-scalar theories of gravity. Particularly, it 
is shown that, under some conditions on the second derivative of the coupling function and on star's 
mass, for a given star there exist two strong-scalar-field solutions as well as the usual weak-field 
one. This last solution happens to be unstable and a star, becoming massive enough to allow for all 
three solutions, evolves to reach one of the strong field configurations. This transition is dynamically 
computed and it appears that the star radiates away the difference in energy between both states (a 
few Mqc^) as gravitational radiation. Since part of the energy (~ 10~^ Mqc^) is injected into 
the star as kinetic energy, the velocity of star's surface can reach up to 10~^c. The waveform of this 
monopolar radiation is shown as well as the oscillations undergone by the star. These oscillations are 
also studied within the slowly-rotating approximation, in order to estimate an order of magnitude 
of the resulting quadrupolar radiation. 

PACS number(s): 04.40.Dg, 04.50. -f-h, 04.80Cc, 04.30. Db 



I. INTRODUCTION 

Tensor-scalar theories have been studied in many works (see e.g. ||l|) as natural generalizations of Einstein's general 
relativity, representing the low-energy limit of superstring theories ( |^ and Q). They describe gravity by the usual 
spin-2 field (g^i/) combined with one or several spin-0 fields {(p). In all this paper, only one scalar field will be 
considered, with no self-coupling (potential). In that case, the most general tensor-scalar theory is given by the 
action: 

S = {WttG,)-' f d^xV^AR* - 2g':^d^^d.^) + ^„[*„„a2((^)5* ], (1.1) 



where all quantities with asterisks are related to the "Einstein metric" g*^: G* is the bare gravitational coupling 
constant, i?* = g*^R*^i, the curvature scalar for this metric and = det((7*^). The term Sm denotes the action of 
matter, represented by the fields ^m, which is coupled to the "Jordan-Fierz" metric g^^, = a?{Lp)g*^j^] all quantities 
with a tilde are related to this metric. That means that all non-gravitational experiments measure this metric, 
although the field equations of the theory are better formulated in the Einstein one. By varying S, one obtains 

R% - \9l.R* = 2d^^d,^ - g;,gpfdp^d,v + ^r;,, (1.2) 

= -47rG,a((^)r*, (1.3) 

where the Einstein-frame stress-energy tensor T^" = 2{—g*)^^^^ SSm/Sg'^^, is related to the physical one, T^^" = 
2{-~g)-^^HSm/6g^^, by 



Ti:^ = aHv)fi:. (1.4) 



The logarithmic derivative of the coupling function a{ip) is a{ip), present in eq. (1.3). It represents the field-dependent 
coupling strength between matter and the scalar field. Hereafter, it is assumed that this coupling strength function 
contains no large dimensionless parameter, and hence the class of coupling functions is well represented by an affine 
one, depending on two parameters (see Q). If t^o denotes the asymptotic value of the scalar field, uq = a{(pQ) and 
Po = da{(po) / dtfio can be chosen]^ to parametrize a{(p): 



^ipo is then redundant with qo and does not represent a parameter for the theory 
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a{ip) = ao + f3o X {(p - ipo). (1.5) 

Brans-Dicke theory is obtained for (Hq = 0, on the opposite, even if ao = (and /3o ^ 0), the scalar field can exhibit 
some non-perturbative effects in neutron stars (strong gravity), when /3q is lower than some critical value (about -5, 
depending on the mass of the star). This phenomenon has been described in Q and ||] as "spontaneous scalarization" , 
by analogy with the spontaneous magnetization of ferromagnets below the Curie temperature. The aim of this paper 
is to determine under which conditions a neutron star shows such spontaneous scalar field and how can the dynamical 
transition be. Resulting gravitational radiation will also be discussed. The organization is as follows: section || 



delimits the parameter space for non-perturbativc effects to appear and sets a maximal value forj3o, section III shows 
computed evolution of unstable stars, during which gravitational radiation is emitted (section IV). Finally, section ^ 
gives some concluding remarks. 



II. CONDITIONS FOR SPONTANEOUS SCALARIZATION 



In all this sect ion, spherical symmetry is assumed (this will not be the case in sec. IV B| ), and the coupling function 
of the form (O) will be considered throughout all the paper. It gives the following expression for the function a{tf): 



a{ip) 



(2.1) 



so that a(</3o) = 1- As stated above, i^jq enters the theory only as a boundary condition for the scalar field, but not as 
a parameter of the theory. 



A. Coordinate and Variable choices 

Following 1^ , space-time is decomposed on spacelike hypersurfaces within the 3 + 1 formalism, with the radial gauge 
and polar slicing. The metric 5*^ takes the form 

ds^ = -N^{r, t)dt^ + A^{r, t)dr^ + r^{de^ + sin^ Odcj)'^) (2.2) 

All coordinates are expressed in the Einstein-frame, and asterisks are omitted. However, "physical" quantities will 
often be written in the Fierz metric and noted with a tilde. Neutron star matter is modeled as a perfect fluid 

where -u^ is the 4-velocity of the fluid, e is the total energy density (including rest mass) in the fluid frame and p is 
the pressure. The description of the fluid is completed by the equation of state, not depending on temperature (cold 
matter) 

e = e(ns), 

with fiB being the baryonic density in the fluid frame. This assumption consists in neglecting the strong and weak 
nuclear interaction processes and assuming matter is at equilibrium for these processes. This is, of course, not valid 
for /3-equilibrium, but as it has been shown by |^], the effects of weak nuclear interactions on the hydrodynamics of a 
neutron star collapse are negligible. Finally, the coordinate velocity is defined as 

dr 

and the Lorentz factor of the fluid as 




The equations for the fields and matter variables are those of [|6| and are solved the same way, by the same numerical 

code. There will be 4 equations of state used in this work: EOSl , a polytrope (p = KfioifiB (^If") ) ^ith 7 = 2.34 and 

K = 0.0195, TOB = 1.66 X 10-2"^ kg and no = 0.1 fm'^; E0S2, the same as above, but with 7 = 2 and = 0.1; E0S3, 
Pandharipande equation of state (realistic, rather soft equation of state; see |^ and for properties) and EOS4, 
Bethe and Johnson equation of state (realistic, rather stiff equation of state; see [|j and [|o| for properties). These 
equations have chosen because they describe different stiffnesses of matter and, on the other hand, their numerical 
behavior allows for good accuracy in the results. 
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B. Maximal Po parameter for Spontaneous Scalarization 



When considering static solutions of eqs. (1.2)-(1.3) for ao = 0, there may exist two types of solutions: one with 
(7* J, being the same as in general relativity, and tp = constant — ipo, the scalar charge of the star is thus null; two 
solutions where \ip — ipo\ ~ 1 at star's center, they are images of each other hy tp —^ —ip + 2ipQ, and have the same 
global characteristics, in those cases, the scalar charge is of the order of the star's mass. 

The first type of solution (which will be noted weak-field solution) always exists, whereas the second one (called 
strong- field solution) requires that the amount of baryons in the star (its baryonic mass Mb) be larger than some 
critical value, depending on the (Jq parameter. This result has been obtained iDy Damour and Esposito-Farese ( [||, 
Fig. 2 and Tab. 1) and Harada ( Tab.l). Using the same equation of state as these two works (EOSl), static 
weak-field solutions were computed for increasing baryonic masses and checked against spontaneous scalarization. 
Since in the ao = case the transition has an infinite slope (Cf. Fig. 1 of |^), it is possible to draw the curve Mg*"** 
as a function of (3o; the results are displayed in Fig. |l|. The curve is parametrized by stars' central density, until it 
reaches its maximal value, corresponding to maximal mass (2.23Mq for general relativity). If this density is further 
increased, neutron stars configurations become hydrodynamically unstable, they may however show spontaneous 
scalarization effects if their mass is above the dotted curve. For hydrodynamically stable neutron stars (being below 
the Mb = 2.23Mq line), strong-field solutions can develop inside the solid curve. This is in accordance with previous 
works ( Q and [pT|) and shows that there exists a maximal (3q parameter for spontaneous scalarization to occur 
(/3o — —4.34). For —5.2 < Po < —4.34, Fig. |l| shows that a star, which is close to its maximal mass, does not exhibit 
spontaneous scalarization effects, whereas for the same /3o less massive stars do. An explanation to this could be that, 
for such stars e — 3p at their center becomes small or negative, therefore the simplified model of static equations of 
1^ shows no more "zero modes" . However, this cannot be the only reason, the (3o parameter should also intervene. 

This study has been extended to other equations of state (E0S2 to 4) and the results are shown in Fig. ^. The 
main result is that the maximal value of (3q for which spontaneous scalarization still occurs is independent from the 
equation of state used. Its value reads 

/J™^'^ = -4.34 ±0.01. (2.3) 

Since the equations of state used cover a very broad range of stiffness and are, for two of them, results of realistic 
scenarios for dense matter, this limit can be considered as a very strong one for compact stars. It can easily be 
interpreted with the results of Harada jl^ : taking an incompressible fluid model, he has shown that the (3o parameter 
below which instabilities of the scalar field develop, is little sensitive to the compactness {R/M) of the star around 



R/M = 4 (Fig. 5 of his work). Section |III| will show that these dynamical instabilities are expression of the existence 
of strong-field solutions. 



C. Spontaneous Scalarization when ao 7^ 

What has been stated in previous section was a result of computations in the case ao = 0. Behavior of neutron 
stars is quite similar in the case ao ^ 0, if we consider this parameter constrained by solar-system experiments (see 
1 and II) to: 

al < 10-^. (2.4) 
The most general conditions for spo ntaneous scalarization to appear have been studied by p^ , using catastrophe 



theory Namely, for ao verifying {2A), there are still two kinds of static solutions: one for which the scalar charge is 
of the order of ao, the second one (containing two different solutions, not equivalent this time) for which the scalar 
charge is of the order of unity. Some models of neutron stars are described in Table || in order to compare models 
with same number of baryons (I.SA/q), with /3o = — 6 for all of them so that spontaneous scalarization is likely to 
occur. One sees that solutions with high scalar charge are energetically more favorable than solutions with no or 
small charge; beside this, spontaneous scalarization mainly affects matter distribution inside the star, but very little 
global variables such as the mass and the radius. One should note that in the case ao ^ 0, when the mass of the 
neutron star increases, the transition to spontaneous scalarization is smoothed and one has to delimit the zone for 
spontaneous scalarization to happen (as in Fig. |^) by looking for the existence of three solutions for a given amount of 
baryons. The study of the number of equilibrium solutions and of their stability has been done, for EOSl, by Harada 
p3[, using catastrophe theory. The same results have been found here, using a static numeric al co de, and have been 



extended to other equations of state, getting thus as a new result the maximal value of Po (2^)i compatible with 
spontaneous scalarization in the ao = case. 



3 



III. DYNAMICAL TRANSITION TO SPONTANEOUS SCALARIZATION STATE 



The conditions for spontaneous scalarization to appear are now well defined and the question is now to know, for a 
fixed theory (i.e. ao) Po and the equation of state) and a given amount of baryons, if the three equilibrium solutions 
are stable and what happens to the unstable ones. The first part of this question has already been answered by Harada 
pTf who did a semi-analytical stability analysis of spherically symmetric neutron stars in tensor-scalar theory; and 
by Harada using the catastrophe theory. He showed the development of unstable modes for weak-field solutions 
which correspond to the possibility for spontaneous scalarization. Here, the dynamical numerical code showed the 
same results. This code is described in more details in Q, the main point being the fact that, thanks to pseudo- 
spectral techniques, it is precise enough to be sensitive to instabilities, and to follow the evolution of the unstable 
star. The stability results have been numerically checked, and the code showed that the weak-scalar-ficld equilibrium 
solutions were unstable when strong-scalar field solutions existed (e.g. solution 1 and 4 of Tabic |). 

This study can be linked with the physical scenario of an accreting neutron star (e.g. in an X-ray binary system), 
which passes the limit of the critical baryonic mass. When it is born, the star may be below this mass and thus 
either in a general relativistic state (if ao = 0) or in a "weak-scalar-field" state. When it passes the critical baryonic 
mass, the star will be in an unstable configuration. This kind of solution was then numerically followed and resulting 
evolution is shown on Figs. || to H, for radius, central density, surface velocity and central scalar field, beginning 
with the unstable solution number 4 of Table |[ One sees that the star undergoes a strong variation of its matter 
distribution, caused by the rapid raise of the scalar field, then behaves like a damped oscillator radiating away its 



kinetic energy through monopolar gravitational radiation (see Sec. |IV A|) and finally settles down to the strong-field 
state (static solution 5 of Table |). One can thus see, in the case ao 7^ 0, the spontaneous scalarization appearing 
dynamically. In the simulations, the star starting form the weak-field configuration, would settle down either to 
the positive strong-field state (solution 5) or to the negative one (solution 6). The final dynamically-evolved fields 
correspond to those obtained from a static code within 1% of error (the code indicates 3%, see 0) and the baryonic 
mass is conserved up to 10"^. 

Such simulations have been done for various masses and coupling function parameters. Results were very similar: 
the weak-field state is always unstable when the two other strong-field states exist; these last solutions are both stable. 
Looking at the energies, one sees that the weak-field solution is a local maximum, whereas spontaneous scalarization 
states are local minima (even if the state having a scalar charge of the opposite sign to ao is energetically less favorable 
than the other one, see table |). The difference in energy is radiated away as monopolar gravitational waves in two 
steps. First, the largest part (about 99%) of the difference in energy between both states is radiated very rapidly 
when the scalar field grows to its new value; then, the rest of the energy is put into the star as kinetic energy which 
is dissipated slowly by the oscillations of the star and the change they induce onto the scalar field. 

This last point is studied more in details. There is no rigorous definition of the kinetic energy of a star in general 
relativity, however if one considers the kinetic energy of a particle of mass m: 

i?L = (F - l)mc2, 

F being its Lorentz factor, then a good choice of Ekin for a star may be given by: 

Ekrn= / 47rF(F-l)gr2dr, (3.1) 

with FeQ replacing mc^ . Ekin is thus computed the same way as in jl^ . On Fig. ^ the evolution of this kinetic 
energy, giving the good Newtonian limit, is plotted and shows well the damped behavior (exponential decay). Such 
neutron stars are really damped oscillators, except for stars with masses close to the critical one (described in Fig. |^) 
which exhibit no oscillations, but only exponential relaxation toward equilibrium state. In any case, the maximal 
velocity reached by star's surface is of the order of 10~^c showing clearly that the appearing strong scalar field has 
an important infiuence on the structure and the hydrodynamics of the star. 

IV. GRAVITATIONAL EMISSION DURING THE TRANSITION 



^the additional Lorentz factor is due to the relativistic contraction of space 
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A. Monopolar component 



Variation of the scalar field during the transition results in an emission of monopolar gravitational waves. Far from 
the star, the scalar field writes: 

^ir,t)=^o + l:F{t-^)+o(^^y (4.1) 

and the interesting "scalar wave" which can be detected by an gravitational wave detector, at a distance d from the 
source, is then related to F through: 

hit) = ^aoFit) (4.2) 

(see eq.(5.6) of with a^ipo) ~ 1). The fimction F{t) (waveform) is plotted in Figs. ^ and || for the rapid variation 
and the damped oscillations respectively. The first step has a characteristic frequency of ^ 200 Hz, whereas the second 
one, lower in amplitude, has ~ 3 kHz. These results depend essentially on the star's mass. From these waveforms, 
one can also deduce the radiated scalar energy, defined as: 

which is plotted in Fig. |l^ as a function of time. The function F{t) is estimated at r = 300 km, see The radiated 
amount corresponds (within errors) to the difference between static cases 4 and 5 of table ^. On Fig. ^ the star 
still radiates some energy after t = 10 ms: it is its kinetic energy of Fig. |^. If one does not consider this oscillations 
(which are small compared to the first raise), the wave represents a transition which is the inverse of that of neutron 
star collapse to a black hole, when the star is strongly scalarized. The scalar field goes from the strong-field value to 
the asymptotic one, (po, whereas in Fig. |^ it goes from (almost) (po to its strong field value. Therefore, the amplitude 
of the wave resulting from a transition of a neutron star to a strong-scalar-field state is very similar to that of a wave 
coming from a neutron star collapse to a black hole. The discussion on this amplitude of the wave made in |^ can be 
applied here to say that the monopolar radiation from such events is unlikely to be observed, with constraints on the 
ao, (3o parameters taken from solar system experiments [Q, if the sources are situated at a larger distance than a few 
100 kpc. This excludes the Virgo cluster galaxies and thus reduces the number of possible sources. 



B. Quadrupolar component 



Prom eq. 4.2 one sees that the ampHtude of the wave interacting with the detector is directly proportional to 
ao- Since cosmological arguments ( Pq] and [0) indicate that the parameter ag should have been driven to by 
cosmological evolution, monopolar waves described in the previous section can be "invisible" to gravitational detectors. 
However, if one considers a slowly-rotating neutron star, which is very close to spherical symmetry, the oscillations 
of its surface (Figs. ^ and ^ will induce a modification of its quadrupolar momentum, the star will emit (usual) 
quadrupolar gravitational radiation, whose detection is not sensitive to ag. This qualitative scenario should also hold 
for rapidly-rotating neutron stars, but the aim of this section is to make an order-of-magnitude crude estimate of the 
emitted quadrupolar wave amplitude. The detection of such quadrupolar waves, with no detection of monopolar ones, 
would give solid constraints on the theory, indicating the possibility for spontaneous scalarization and constraining 
ao- 

Therefore, the procedure is very simplified: at each time-step of the dynamical integration the star is supposed to 
be at equilibrium. As a slowly-rotating star, it can be described by the perturbation equations of Hartle p7[ | and 
Hartle and Thorne , which are second order accurate in the angular velocity. Their work starts from a spherically 
symmetric static star, described in the same gauge as in this work, which is perturbed by the rotation. The resulting 
metric involves several additional functions, proportional to angular velocity or to its square: 



-N^[l + 2{ho+h2P2)W 



1 + -A'^{ma+m2P2) 



^[1 + 2(v2 - h2)P2][d9^ + sin^ 0{d(l> - ujdtf] + 0{n^), 



(4.4) 



with il bei ng star's angular velocity, oj{r) the "angular velocity of the local inertial frame", A and TV are defined 
in cq. (2.2), P2 = P2{cos9) is the second Legendre polynomial and /iq, /i2, "Iq, m2, W2 are perturbation functions 
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proportional to il^. From the asymptotic behavior of these last functions, one can deduce the quadrupolar momentum. 
Since all the quantities are computed at each time-step, from the matter distribution and fields, one also knows the 
quadrupolar momentum Q of the star at each time-step. One then uses the usual quadrupole formula to get the 
amplitude of the wave at a distance d: 

2G. d^Qjt) 



The value obtained for Q from the simulation described in sec. Ill is 



^^lO^^x/^kg.m^s-^ 
with / being the rotation frequency of the star, expressed in Hz. This yields 

Ktt ^ 10"^^ X (4.6) 

So, even if the highest value for the rotation frequency is considered (2 kHz, see [^), if this wave is to be detected 
by the interferometers currently under construction (whose detection level is about 10"^'^ at these frequencies), the 
source has to be situated closer than 10 pc, that is very close in our Galaxy. 



V. CONCLUSIONS 



Spontaneous scalarization effects appear in a restricted part of the {Pq,Mb) parameter space, depending on the 
equation of state used for neutron star matter. On the contrary, the maximal Po parameter which allows for spon- 
taneous scalarization is quasi-independent of the equation of state, its value being —4.34. When these types of 
equilibrium solutions exist, they are energetically more favorable than the weak-field ones, which are then unstable 
equilibria. A star which is in the weak-field regime (when its mass is to low to get "scalarized" ) can increase the 
mass by accretion and thus become unstable. Then it undergoes evolution to the strong-field state radiating the 
difference of energy as monopolar gravitational wave; it also behaves like a damped oscillator, since a small fraction 
of the energy (~ 10~^Mqc^) is put as kinetic oscillating energy into the star. The transitions studied in this work 
are not likely to be observed through their gravitational radiation, excepted if they are located close to us (which 
severely decreases the number of possible sources). However, the kinetic energy put into the star has an important 
influence on the hydrodynamics and its effect must not be neglected when considering the supernovae collapse toward 
a neutron star: when the proto-neutron star gets its final compactness the scalar field can develop and influence the 
ejection of the envelop. The future study of supernova collapse and bounce, within the framework of tensor-scalar 
theory could provide us with monopolar gravitational signals. These latter can give constraints on the tensor-scalar 
theory space parameter, even if they are not detected, since in the case of a supernova, electromagnetic or neutrino 
signals are detected. 
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FIG. 1. Zone of spontaneous scalarization of general-relativistic neutron star solutions in the baryonic mass (Mb) ~ l3o 
plane, for EOSl. The zone lies inside solid lines, the horizontal line at 2.23Mq represents the maximal mass for neutron 
stars in general relativity. Since the curve is parametrized by star s' ce ntral density {ns), it has been continued for unstable 
configurations (thin dotted line) for which dMs/dnB < 0. See Sec. [IB for more explanations 
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FIG. 2. Same as Fig. but for EOS2, EOS3 and EOS4 (see Sec. II A) 
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FIG. 3. Evolution of star's radius for unstable initial equilibrium solution number 4 (Cf. Table 
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FIG. 5. Evolution of star's surface velocity for unstable initial equilibrium solution number 4 (Cf. Table 
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FIG. 6. Evolution of star's central scalar field value for unstable initial equilibrium solution number 4 (Cf. Table 
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FIG. 7. Evolution of star's kinetic energy defined as J'^ ^ 47rr(r — l)er^dr, where e is the total energy density of fluid, for 
unstable initial equilibrium solution number 4 (Cf. Table 'ffi . 
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FIG. 8. Scalar waveform resulting from the evolution of unstable initial equilibrium solution number 4 (Cf. Table 
corresponding to the raise of the scalar field. 
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FIG. 9. Scalar waveform resulting from the evolution of unstable initial equilibrium solution number 4 (Cf. Table 
corresponding to the damping of star's oscillations. 



12 



I 

O 
X 

" ' tr> 
^ I 
o o 

© ^ 



ro 

o I 

CO C_J 



X 

CNJ 



o 



6 8 10 12 14 

time [ms] 

FIG. 10. Evolution of star's radiated energy defined in eq.( [l.3| ), for unstable initial equilibrium solution number 4 (Cf. 
Table §) . 
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TABLE I. Computed characteristics of different neutron stars of 1.5Mq showing spontaneo us sc alarization effects (solutions 
2,3,5 and 6) or without (solutions 1 and 4). The equation of state (EOSl) is described in sec. [I A, ao and /3o are the coupling 
function parameters (1.5), Rstar denotes star's radius, hB(r — 0) is the central baryon density (in units of nuclear density. 



1 nnuc ~ 10** m ^), Mq is the (7*,^-frame ADM mass, Mb the baryonic one and lo the scalar charge. The code used for these 
results is a static one. 



Solution 




Po 


^star 


ns(r = 0) 


Ma 


Mb 










[km] 




[Me] 


[Mq] 


[Mq] 


1 





-6 


13.2 


3.9643 


1.37803 


1.50009 





2 





-6 


13 


5.3212 


1.37322 


1.50008 


0.781 


3 





-6 


13 


5.3212 


1.37322 


1.50008 


-0.781 


4 


10"^ 


-6 


13.2 


3.9742 


1.37807 


1.50007 


-5.91 X 10-2 


5 


10"^ 


-6 


13 


5.3674 


1.3719 


1.50008 


0.803 


6 


10"^ 


-6 


13 


5.2669 


1.37452 


1.50008 


-0.757 
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